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Abstract 

The real and imaginary part of the vacuum polarisation function n(q 2 ) induced by 
a massive quark is calculated in perturbative QCD up to order a 2 s . The method is 
described and the results are presented. This extends the calculation by Kallen and 
Sabry from two to three loops. 
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1 Introduction 



The measurement of the total cross section for electron positron annihilation into hadrons 
allows for a unique test of pertubative QCD. The decay rate T(Z — *> hadrons) provides one 
of the most precise determinations of the strong coupling constant a s . In the high energy 
limit the quark masses can often be neglected. In this approximation QCD corrections to 
R = cr(e + e~ — > hadrons)/a(e + e~ — > have been calculated up to order [I], 0. 

For precision measurements the dominant mass corrections must be included through an 
expansion in m 2 /s. Terms up to order a z s m 2 / s || and a^m^/s 2 M are available at present, 
providing an acceptable approximation from the high energy region down to intermediate 
energy values. For a number of measurements, however, the information on the complete 
mass dependence is desirable. This includes charm and bottom meson production above 
the resonance region, say above 4.5 GeV and 12 GeV, respectively, and, of course, top 
quark production at a future electron positron collider. 

To order a s this calculation was performed by Kiillen and Sabry in the context of QED 
a long time ago |JJ . With measurements of ever increasing precision, predictions in order 
a 2 are needed for a reliable comparison between theory and experiment. Furthermore, 
when one tries to apply the 0(a) result to QCD, with its running coupling constant, the 
choice of scale becomes important. In fact, the distinction between the two intrinsically 
different scales, the relative momentum versus the center of mass energy, is crucial for a 
stable numerical prediction. This information can be obtained from a full calculation to 
order a 2 only. Such a calculation then allows to predict the cross section in the complete 
energy region where perturbative QCD can be applied — from close to threshold up to 
high energies. 

In this paper results for the cross section are presented in order a 2 . They are obtained 
from the vacuum polarisation II (q 2 ) which is calculated up to three loops. The imaginary 
part of the "fermionic contribution" - - derived from diagrams with a massless quark 
loop inserted in the gluon propagator — has been calculated in ||. All integrals could 
be performed to the end and the result was expressed in terms of polylogarithms. In this 
paper the calculation is extended to the full set of diagrams relevant for QCD. Instead of 
trying to perform the integrals analytically, we use the large q 2 behaviour of n(g 2 ) up to 
terms of order m 2 /q 2 and calculate its Taylor series around q 2 = up to terms of order 
q 8 . The leading and next-to- leading singularity is deduced from the known behaviour 
of the nonrelativistic Green function and the two-loop QCD potential. Altogether eight 
constraints on II(g 2 ) are thus available, four from q 2 = 0, two from q 2 — > — oo and two 
from the threshold. 

The contributions ~ C F , ~ CaCf and ~ CfTtii have to be treated separately since 
they differ significantly in their singularity structure. For each of the three functions 
an interpolation is constructed which incorporates all data and is based on conformal 
mapping and Pade approximation suggested in J7|, ||, ||, [10| . Since the result for CpTrii 
is available in closed form the approximation method can be tested and shown to give 
excellent result for this case. Reliable predictions for R to order a 2 and arbitrary m 2 / s 
are thus available. 

In this paper only results without renormalization group improvement and resum- 
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mation of the Coulomb singularities from higher orders are presented. Resummation of 
leading higher order terms, phenomenological applications and a more detailed discussion 
of our methods will be presented elsewhere. 



2 The Calculation 

A large number of ingredients are needed in the calculation of the three-loop vacuum po- 
larisation function: The high energy behaviour of II(g 2 ), its Taylor series at q 2 = and its 
singularities at threshold which are related to the QCD potential. The decomposition of II 
according to its colour structure and the separation of gluonic and fermionic contributions 
is crucial to display the striking differences in the threshold region. The approximation 
method is based on conformal mapping and Pade approximation. In the present approach 
both real and imaginary parts of the vacuum polarization II (q 2 ) are calculated and its 
analyticity properties are exploited heavily. 

The physical observable R(s) is related to IT(g 2 ) by 

R(s) = 127rimn(g 2 = s + ie). (1) 

It is convenient to define 

II(g 2 ) = n(°)(g 2 ) + ^!V)(g 2 ) + ( ^-l) 2 T[( 2 \q 2 ) + ... (2) 

7r \ n J 

with the MS coupling a s (n 2 ) defined in the conventional way. 

To describe the singularity structure of II in the region close to threshold the pertur- 
bative QCD potential fill 



W? 2 ) = -47rC F ^4p, (3) 

q 

will become important. The CaCf and CfTui contributions have been displayed sepa- 
rately, with Cf = 4/3, Ca = 3 and T = 1/2. The number of fermions is denoted by Nf 
and has to be distinguished from the number of light (= massless) fermions n\ = Nf — 1. 
To transform the results from QED to QCD the proper group theoretical coefficients 
Cf = 1, Ca = and T = 1 have to be used. 

In this paper we are only concerned with contributions to II(g 2 ) and R(s) which 
originate from diagrams where the electromagnetic current couples to the massive quark. 
In order a s and a 2 all these amplitudes are proportional to Q 2 , the square of the charge 
of the massive quark. Diagrams where the electromagnetic current couples to a massless 
quark and the massive quark is produced through a virtual gluon have been calculated in 
|T2|1 and will not be discussed here. 

For the a 2 calculation the following steps have been performed: 
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Decomposition according to the colour structure 

The contributions from diagrams with n\ light or one massive internal fermion loop 
will be denoted by CpTntH) and CpTYi F with the group theoretical coefficients factored 
out. Purely gluonic corrections are proportional to C\ or CaCf- The former are the only 
contributions in an abelian theory, the latter are characteristic for the nonabelian aspects 
of QCD. It will be important in the subsequent discussion to treat these two classes 
separately, since they exhibit qualitatively different behaviour close to threshold. The 
following decomposition of II (q 2 ) (and similarly for R(s)) is therefore adopted throughout 
the paper 

n = n( > + ^^n« (5) 

C A C F U^ A + CpTnfr® + C F TUP] . (6) 

All steps described below have been performed seperately for the first three contributions 
to n< 2 >. In fact, new information is only obtained for and since ImIIp' ) is known 
analytically already. The amplitude with a massive internal fermion exhibits a two particle 
cut with threshold at 2m which has been calculated analytically || . The contribution from 
a four particle cut with threshold at 4m is given in terms of a two dimensional integral |J 
which can be solved easily numerically. The II)/ term will not be treated in this paper. 

Analysis of the high q 2 behaviour 

The high energy behaviour of II provides important constraints on the complete an- 
swer. In the limit of small m 2 /q 2 the constant term and the one proportional to m 2 /q 2 
(modulated by powers of ln// 2 /g 2 ) have been calculated a long time ago [p"3[] . For the 
imaginary part even the m 4 /g 4 terms are available [|J. This provides an important test of 
the numerical results presented below. 

Threshold behaviour 

General arguments based on the influence of Coulomb exchange close to threshold, 
combined with the information on the perturbative QCD potential and the running of 
a s dictate the singularities and the structure of the leading cuts close to threshold, that 



is for small v — y 1 — Am 2 / s. The Cp term is directly related to the QED result with 
internal photon lines only. The leading 1/v singularity and the constant term of Ra can 
be predicted from the nonrelativistic Greens function for the Coulomb potential and the 
0(a s ) calculation. The next-to- leading term is determined by the combination of one loop 



results again with the Coulomb singularities [14, 15, 10]. One finds 




Rf = 3^-3^ + ... . (7) 



The contributions ~ CaCf and ~ CpTrii can be treated in parallel. The leading 
CaC f and CpTni term in R is proportional to lnv 2 and is responsible for the evolution 
of the coupling constant close to threshold. Also the constant term can be predicted by 
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the observation, that the leading term in order a s is induced by the potential. The 0{a s ) 
result 

* - »(r +0 'xr + -) < 8 > 

is employed to predict the logarithmic and constant CfCa and CpTni terms of 0(a 2 ) by 
replacing a s by av(4p 2 = v 2 s) as given in Eq.([|). This implies the following threshold 
behaviour: 

(2) ,/31 11 , 11 u 2 \ 

R?> = 3^(-l + il„^ + iln^ + ...). (10) 

This ansatz can be verified for the CfTui term since in this case the result is known in 
analytical form |J . Extending the ansatz from the behaviour of the imaginary part close 
to the branching point into the complex plane allows to predict the leading term of IT(g 2 ) 
~ \nv and ~ lnt> 2 . 

Behaviour at q 2 = 

Important information is contained in the Taylor series of IT(g 2 ) around zero. The 
calculation of the first four nontrivial terms is based on the evaluation of three-loop 
tadpole integrals with the help of the algebraic program MATAD written in FORM [plj 
which performs the traces, calculates the derivatives with respect to the external momenta. 
It reduces the large number of different integrals to one master integral and a few simple 
ones through an elaborate set of recursion relations based on the integration-by-parts 



method fL7| , |lq| . Though this master integral, calculated in jL8|, appears in the result 
for single diagrams, in the final expession it cancels. The following results are obtained 
(excluding the H F term): 

8687 w , 32 w , , 22781 



F V 864 sv ; 5 SV ; 1728 SV V 



... , , 142 328, /r !(» , 

+ CpTnt In At 1(2) 

1 243 243 m 2 15 s V 



/ 223404289 24 *92 2 I 4857587 ^.l 

° 2 ~ ° F [~ 1866240 + Y C(2) "^5" C(2)ln2 + ^608T C(3) 

„ „ ( 1030213543 4939, u 2 32 _ . 96 , , « , 723515 > , ' 

+ C A C F + In At ( 2 + — ( 2 In 2 + ( 3 

\ 93312000 2025 m 2 35 sv 7 35 sv 7 55296 sv , 

„ m ( 40703 1796, u 2 32^ A 

+ CfTui In -A C 2 , 

V 60750 2025 m 2 35 SV ; / ' 
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= /_ 885937890461 64 _ 512 33067024499 

F \ 1161216000 21 sv ; 105 w 51609600 w 

„ „ ( 95905830011197 2749076, u 2 256 A/ , 256 A/ 
+ C A C F \ + In At C(2) + C 2 ln2 

V 1706987520000 1488375 m 2 315 s v ; 105 s v ; 

5164056461 _ A „ m ( 9703588 999664 , fi 2 256 A/ A 

+ C(3))+C F Tni\ In At C(2) , 

103219200 sv V V 17364375 1488375 m 2 315 sv 7 ' 

2 / 269240669884818833 640 1024 1507351507033 

C4 — Cp 1 C(2) — C(2) ln2 H C(3) 

F V 61451550720000 231 w 231 w 412876800 w 

„ „ ( 36675392331131681 571846, a 2 512 >/ , 512 >/ 
+ CaCf + In -A C(2) + C 2 ln2 

V 158018273280000 382725 m 2 693 s v ; 231 s v ' 

1455887207647 >/ A ^ m / 54924808 207944, u 2 512 >/ A 

+ C(3) ) + CpTni ( In -A C(2 . 

7431782400 sv 7 \ 120558375 382725 m 2 693 s v 7 

The terms in Ci, C2 and C3 are in agreement with [10|, the remaining ones are new. 



(Note that the C F T contribution calculated in JTD| were obtained with massive internal 
fermions.) 

Conformal mapping and Pade approximation 

The vacuum polarisation function II^ 2 ^ is analytic in the complex plane cut from q 2 = 
Am 2 to +00. The Taylor series in q 2 is thus convergent in the domain \q 2 \ < Am 2 only. 
The conformal mapping which corresponds to the variable transformation 



Izilz q2/Am2 g 2 _ (12) 

1 + Jl-g 2 /4m 2 ' 4m 2 (1 + w) 2 ' 



transforms the cut complex q 2 plane into the interior of the unit circle. The special points 
q 2 = 0, 4m 2 , —00 correspond to uj = 0, 1, —1, respectively. 

The upper (lower) part of the cut is mapped onto the upper (lower) perimeter of 
the circle. The Taylor series in uj thus converges in the interior of the unit circle. To 
obtain predictions for n(g 2 ) at the boundary it has been suggested 0, || to use the Pade 
approximation which converges towards n(g 2 ) even on the perimeter. To improve the 
accuracy the singular threshold behaviour and the large q 2 behaviour is incorporated into 
simple analytical functions which are removed from before the Pade approximation 
is performed. The quality of this procedure can be tested by comparing the prediction 
with the known result for Imnp' ) . 

The logarithmic singularities at threshold and large q 2 are removed by subtraction, 
the 1/v singularity, which is present for the Cp terms only, by multiplication with v as 
suggested in fLOf . The imaginary part of the remainder which is actually approximated 
by the Pade method is thus smooth in the entire circle, numerically small and vanishes 
at uj — 1 and uj = —1. 
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3 Results 



After performing the Pade approximation for the smooth remainder with u as natural 
variable, the transformation (12) is inverted and the full vacuum polarisation function 
reconstructed by reintroducing the threshold and high energy terms. This procedure pro- 
vides real and imaginary parts of IT 2 ). Subsequently only the absorbtive part of II^ 2 ^ 
(multiplied by 12n) will be presented. 

In Figure [I] the complete results are shown for fi 2 = m 2 with R ( 2\ R ( na and Rf ] 
displayed separately. The solid curves are based on the [4/2], [3/2] and [3/2] Pade ap- 
proximants for A, NA and I, respectively. The threshold and high energy behaviour is 

(2) 

given by the dashed curves. The exact analytical result which is known for the R\ con- 
tribution only differs from the approximation curve in Figure Jl] by less than the thickness 
of the line. The quality of the approximation for and R^ is confirmed by the com- 
parison of the high energy behaviour of the approximation with the known asymptotic 
behaviour (Figure 0). The quadratic approximation (dash-dotted line) is incorporated 
into R^ by construction, the quartic approximation shown (dashed line) is known from 
HI, but is evidently very well recovered by the method presented here. 

Different Pade approximations of the same degree ([4/2] or [2/4]) and approximants 
with a reduced number of parameters give rise to practically indentical predictions, which 
could hardly be distinguished in Figures [I] and 0. Minor variations are observed close 
to threshold, after subtracting the singular and constant parts. The remainder 6R for 
three different Pade approximants is shown in Figure |3|. The stability of the different 
approximations and their smooth behaviour close to the point v = 1 can be considered as 
additional confirmation of our ansatz in Eqs.([7y£||i~0|). The spread of the different curves 

(2) 

can be used to estimate the quality of the approximation. For R) the approximation 
can be compared to the exact result (dotted line) and the nearly perfect agreement gives 
additional support to the approach presented in this work. 

To summarize: Real and imaginary part of the vacuum polarisation function H(q 2 ) 
from a massive quark have been calculated up to three loops for QCD and QED. This 
result extends the classic calculations of Kallen and Sabri || to next-to-leading order. 
The imaginary part can be used to predict the cross section for production of massive 
quarks for arbitrary m 2 / s, wherever perturbative QCD can be justified — from above the 
quarkonium resonance region up to high energies. 



Acknowledgments 

We would like to thank A.H. Hoang and T. Teubner for many interesting discussions. Their 
programs with the analytical results for R) were crucial for our tests of the approximation 
methods. We are grateful to D. Broadhurst for helpful discussions. J.H.K. would like 
to thank S. Brodsky for discussions of the V-scheme and the threshold behaviour and 
K.G.Ch. thanks A. A. Pivovarov for careful reading the manuscript. 



6 



References 

[1] K.G. Chetyrkin, A.L. Kataev, F.V. Tkachov, Phys. Lett. B85 (1979) 277; 
M. Dine, J. Sapirstein, Phys. Rev. Lett. 43 (1979) 668; 
W. Celmaster, R.J. Gonsalves, Phys. Rev. Lett. 44 (1980) 560. 

[2] S.G. Gorishny, A.L. Kataev, S.A. Larin, Phys. Lett. B259 (1991) 144; 

L.R. Surguladze, M.A. Samuel, Phys. Rev. Lett. 66 (1991) 560; erratum ibid, 2416. 

[3] K.G. Chetyrkin, J.H. Kiihn, Phys. Lett. B248 (1990) 359. 

[4] K.G. Chetyrkin, J.H. Kiihn, Nucl. Phys. B432 (1994) 337. 

[5] G. Kallen and A. Sabry, K. Dan. Vidensk. Selsk. Mat.-Fys. Medd. 29 (1995) No. 17, 
see also J. Schwinger, Particles, Sources and Fields, Vol.11, (Addison- Wesley, New 
York, 1973). 

[6] A.H. Hoang, J.H. Kiihn, T. Teubner, Nucl.Phys. B452 (1995) 173. 

[7] J. Fleischer and O.V. Tarasov, Z. Phys. C64 (1994) 413. 

[8] D.J. Broadhurst, J. Fleischer and O.V. Tarasov, Z. Phys. C60 (1993) 287. 

[9] D.J. Broadhurst et. al., Phys. Lett. B329 (1994) 103. 

[10] PA. Baikov and D.J. Broadhurst, Report Nos. OUT-4102-54, INP-95-13/377 and 
|hcp-ph/9504398| , to appear in New Computing Techniques in Physics Research IV, 
(World Scientific, in press). 

[11] W. Fischler, Nucl. Phys. B129 (1977) 157. 
A. Billoire, Phys. Lett. B92 (1980) 343. 

[12] A.H. Hoang, M. Jezabek, J.H. Kiihn, T. Teubner, Phys. Lett. B338 (1994) 330. 

[13] S.G. Gorishny, A.L. Kataev and S.A. Larin, Nuovo Cimento 92 A (1986) 119. 

[14] R. Barbieri, R. Gatto, R. Kogerler and Z. Kunszt, Phys. Lett. B57 (1975) 455. 

[15] B. H. Smith and M.B. Voloshin, Phys. Lett. B324 (1994) 117. 

[16] J. A.M. Vermaseren, Symbolic Manipulation with FORM, (Computer Algebra Nether- 
lands, Amsterdam, 1991). 

[17] F.V. Tkachov, Phys. Lett. B100 (1981) 65. 

K.G. Chetyrkin and F.V. Tkachov, Nucl. Phys. B192 (1981) 159. 

[18] D.J. Broadhurst, Z. Phys. C54 (1992) 54. 



7 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

v 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

v 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

v 



Figure 1: Complete results plotted against v = \Jl — 4m 2 / 's. The high energy approxi- 
mation includes the m 4 /s 2 term. 
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Figure 2: High energy region. The complete results (full line) are compared to the high 
energy approximations including the m 2 / s (dash-dotted) and the m A / s 2 (dashed) terms. 
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Figure 3: Threshold behaviour of the remainder SR for three different Pade approximants. 
(The singular and constant parts around threshold are subtracted.) 
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